Data compression of measurements of peculiar velocities of Supernovae la 
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We study the compression of information present in the correlated perturbations to the luminosity 
distance in the low-redshift (z < 0.1) supernovae la due to peculiar velocities of these supernovae. 
We demonstrate that the naive compression into angular velocity power spectrum does not work 
efficiently, due to thickness of the spherical shell over which the supernovae are measured. Instead, we 
show that measurements can be compressed into measurements of f 2 P(k), where / is the logarithmic 
rate of growth of linear perturbations and P(k) is their power spectrum. We develop an optimal 
quadratic estimator and show that it recovers all information for ACDM models for surveys of 
N ~ 10, 000 or more supernovae. We explicitly demonstrate robustness with respect to the assumed 
fiducial model and the number of power spectrum bins. Using mock catalogues of SNe la we estimate 
that future low redshift surveys will be able to probe erg to 6% accuracy with 10, 000 SNe la. 
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I. INTRODUCTION 

Since the earliest studies of supernovae, it has been 
suggested that type la (SNe la) might be used as stan- 
dard candles for cosmological measurements. In the sub- 
sequent years many studies using SNe la revealed that 
the expansion of the universe is accelerating p] [2] . From 
SNe la observational data one could measure cosmolog- 
ical parameters describing the homogeneous expansion 
of the universe through measurements of the luminosity 
distance such as matter and dark energy densities and 
equation of state parameters [5HS] • 

But density inhomogeneities cause additional scatter 
to SNe Hubble flow [10-23 . In the recent years it has 
become apparent that the correlations between these pe- 
culiar velocities of the supernovae are significant at low 
rcdshifts (z < 0.1) [24l [25] and thus present an opportu- 
nity to measure cosmological parameters that affect the 
growth of perturbations in the Universe such as ampli- 
tude and shape of the matter power spectrum. 

But extracting any additional information from the 
correlations induced by peculiar-velocities comes at an 
expense of computing N sn x N sn correlation matrix and 
its inverse on every likelihood evaluation. To complicate 
matters even further, each matrix element consists of an 
oscillatory integrals thus adding to the cost of the over- 
all computing time. In this paper we asses whether a 
data compression method can be found that will provide 
both efficiency and simplified form of the cosmological 
information carried by the peculiar velocity correlation 
matrix. This will also provide a physical insight into 
what is that the peculiar velocities of supernovae la are 
measuring. 

Through the paper we assume the correlated part of 
peculiar velocities is described by the linear theory and 
that small scale virial velocities can be described by sin- 
gle parameter describing the velocity dispersion. While 
this is probably not a good approximation for realistic 
surveys of thousands of supernovae (after all, supernovae 
don't measure velocities at random positions and the real 



Universe is not linear), it nevertheless provides a good 
simple framework for studying an efficient information 
compression. 

The paper is structured as follows. In Section [IT] we 
discuss the background theory required to describe corre- 
lations in the luminosity distances to nearby supernovae 
sourced by correlated large-scale velocities. In Section [TTT| 
we discuss the angular power spectrum of the projected 
velocities as a candidate for data compression. In Section 
IV we discuss the power spectrum multiplied by a factor 
describing the growth of fluctuations at mean redshift as 
our proposed method for data compression. We conclude 
in Section IVl 

Throughout the paper we assume that unless noted 
otherwise, the variables have their usual meaning. 



II. CORRELATIONS OF PECULIAR 
VELOCITIES 

The luminosity distance to a distant SN at redshift 
z is defined by 
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where F is the measured flux and L the intrinsic lumi- 
nosity of the supernovae. In the Friedmann-Robertson- 
Walker (FEW) metric, the luminosity distance is related 
to the comoving distance to an object at redshift z as 
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where k is the geometrical curvature of the universe and 
we have introduced x{ z ) as the comoving distance to a 
distant object at z as 
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Astronomers prefer the flux relation (Eq. ([I 
in terms of the magnitudes as 
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where m and M stand for the apparent and absolute 
magnitude respectively. 

The perturbed FRW metric leads to the effect of pecu- 
liar velocities and those in turn lead to the perturbations 
in luminosity distance given by 
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where d™ e and stand for measured luminosity dis- 
tance and theoretical prediction for unperturbed space- 
time, given by Eq. (|2). Velocities v and v Q are pecu- 
liar velocities of the source and the observer respectively. 
Projection of peculiar velocities along the line of sight (x) 
is the only component we can measure. Redshift z on the 
right-hand side of the equation stands for the observed 
redshift. 

Since the peculiar velocities result from some ini- 
tial Gaussian matter perturbations, the peculiar veloc- 
ity measurements are drawn from a distribution with 
zero mean and nonzero variance. The later can be writ- 
ten in the form of correlation function as £(xi,X2) = 
((v(xi) • xi)(v(x2) • X2)*) for two SN at comoving po- 
sition Xi 2- The correlation function of projected pecu- 
liar velocities has been computed in a number of studies 
[2"S"ll27| . The result, assuming linear perturbation theory, 
can be written as 
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where x 12 = X! — x 2 , x = |x 12 |, cos 6^ = Xi • x 12 and 
cos 62 = x 2 • xi 2 . The projections of the correlation func- 
tion are given by [2FJ 
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where D(z) is the normalized growth function (D(z) = 
D(z)/D(z = 0)) and derivatives are with respect to con- 
formal time (' = d/dj], a(t)drj = cdt). P(k) is the mat- 
ter power spectrum and integration kernels are given by 

Correlations in the peculiar velocity lead to correla- 
tions in the luminosity distance fluctuations which can 
be written as [23 
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The total correlation matrix is a sum of peculiar velocity 
correlation matrix and uncorrelated scatter of and can 
be written as [251 



Cij = C v (i,j) + S tJ af. 
The diagonal parts of the matrix given by |25j 



(9) 



In 10 



.)+ 1 



(1 + *)' 
H{z)d L 



(10) 



where a mi stands for observational errors on apparent 
magnitudes, which are a property of the dataset and 
vary from supernova to supernova. The remaining two 
parameters are a m , which is intrinsic magnitude scatter 
describing deviations of supernova luminosities from the 
perfect standard candles, while a v is the small scale ve- 
locity dispersion due to uncorrelated small scale peculiar 
velocities. 



III. PROJECTED VELOCITY ANGULAR 
POWER SPECTRUM 

Measurements of the SNe are most commonly trans- 
formed into luminosity distances or, in case of pecu- 
liar velocities, into luminosity distance fluctuations. But 
there is no apparent reason why not to work with the 
projected peculiar velocity field instead [22]. Since we 
are interested in the velocity perturbations along the line 
of sight, the projected peculiar velocity field seems more 
natural as well. Those velocity perturbations can be un- 
derstood as perturbations relative to the expansion of the 
universe in physical coordinates, or as perturbations with 
respect to the comoving grid. 

Since the peculiar velocity field is a scalar function 
on a sphere, it is natural to consider its angular power 
spectrum. Using expansion over spherical harmonics we 
get 
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where we kept the dependence of the expansion on red- 
shift z, i.e. we consider infinitely thin spherical shells 
at redshift z. Yg m are spherical harmonic functions that 
satisfy orthonormal relation 
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Using the assumption of Gaussian initial perturbations 
we can relate angular auto and cross power spectra to 
the expansion coefficients (a^ m ) as 
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FIG. 1: The luminosity distance fluctutations angular 
power spectrum plotted as £{£ + l)De/2ir vs £ for three 
different redshifts: z = 0.05 (solid), z — 0.1 (dashed) 
and z = 0.01 (dot-dashed). The angular power 
spectrum has dimensionless units. 



FIG. 2: The cross correlation coefficient r(z) (Eq. 16) 
with respect to reference redshift z m = 0.016528. The 
coefficient is plotted for four low multipoles: £ = 2 
(solid), £ = 3 (dashed), £ = 4 (dot-dashed) and £ = 5 
(dotted). 



In the linear theory, these are given by 
C i (z,z')=D , (z)D'{z')-- 



(14) 



for two SNe at the redshifts z and z' 7 where D is the 
growth factor, ji are spherical Bessel functions, x ' = 
x(z( )) and P(k) matter power spectrum. The detailed 
derivation can be found in the Appendix. 

The angular power spectrum of luminosity distance 
fluctutations are exactly the same as in Eq. ( 14 ) save 
for the prefactors and can be written as 
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Both angular power spectra have dimensionless units. 
This is due to the fact that velocity in our all calculations 
is in the units of the speed of light. 

Figure ([I I shows the auto-power spectrum Di(z, z) for 
different values of redshift. As expected and also shown 
in the literature [22], numerical simulations are in good 
agreement with above result (Eq. (14|) only for small £ 
(£ < 100), where linear regime is still valid. 

Although the formalism is appealing, is the angular 
power spectrum actually a useful compression method 
for the low-redshift supernovae? If we want to describe 
all peculiar velocities of supernovae at low redshifts using 
a single power spectrum, the slices at different redshifts 
would essentially need to be describing the same velocity 



field. In order to check this assumption, we look at the 
cross-correlation coefficient given by 



r(z,z') = 



C}(z,z>) 



C e (z,z)C e {z',z>) 



(16) 



In Figure ^ we plot r(z, z') while fixing z' — 0.016, 
a median redshift of a typical low-z SNe dataset used in 
|17j . We see that the 'correlation length' of velocities in 
redshift is only about Sz = 0.01 and hence to describe all 
correlations between z — 0.0 and z = 0.1, we would need 
about 10 angular power spectra, together with a large 
covariance matrix that will describe not only the errors 
and their covariances at one redshift, but also those at 
neighbouring redshifts. This makes this approach clearly 
suboptimal and hence we turn the direct estimation of 
the three-dimensional power spectrum. 



IV. EFFECTIVE POWER SPECTRUM 
MEASUREMENTS 

In any data compression technique, we try to put con- 
straints on quantities that are as close to the data as 
possible and as independent of the underlying theoreti- 
cal assumptions, while at the same time still capturing 
most of the information in the full dataset. Although 
the covariance matrix of perturbations to the luminosity 
distance depends on a number of cosmological dependent 
pre-factors, in addition to the underlying power spectrum 
P(k), we show that, in the limit of supernovae being 
at sufficiently low redshifts, the data effectively measure 
P(k)f 2 , where / is the logarithmic growth rate, regard- 
less of which fiducial cosmology one adopts. 

Since the P(k) is a continuous function of the wave 
vector k we will approximate it with a stepwise function 
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such that P(k) — const, for k a < k < k a+ \. With this 
in mind we can rewrite the parallel in perpendicular pro- 
jections of correlation function of two SNe (i and j) as 
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where x = |xy|, = Xj-Xj and Xj = x{ z i)*i (Eq. ((3)), 
D(z) is normalized growth function and j± are spherical 
Bessel functions of the first kind. We have also defined 
integral sinus as Si(x) = J* m j^dz. Eqs. (17) and (18) 
are valid only when x 7^ 0. Because x is the norm of 



the difference of two vectors to SN we will surely have 
examples where x = (peculiar velocity auto-correlation 
functions). When taking x — > the above expressions 
simplify to 

£||,x([*i,*2)) = &{zi)&(zj)?^- (k 2 - h) . (19) 

Parameters we want to constrain from the SNe la data 
are values of matter power spectrum P a , and two uncor- 
rected errors a m and a v and so N p = Nb + 2. We have 
assumed that the absolute magnitude offset will be for all 
our purposes completely constrained from all supernovae 
data (including higher-redshift ones). 



A. Mock catalogues 

To test our method we created several mock catalogues 
of SNe redshifts and their positions on the sky. For 
those catalogues we calculated their synthetic observa- 
tional data in the form of luminosity distance fluctua- 
tions. We chose to model are synthetic data with fidu- 
cial cosmological parameters (O m = 0.24, ft, = 0.7, w — 
— 1,0a = 1 — Qm)- We computed the covariance matrix 
give by Eq. Q using reference linear matter power spec- 
trum P re f (k) computed using CAMB 28J. We computed 
eigenvalues A; and eigenvectors v\ { of covariance matrix 
and than, for each A;, draw a random number from a 
Gaussian distribution with a = a/A7 and add this num- 
ber to the initially zero data vector in the direction of 
the V\ t . This procedure gave us at the end a data vec- 
tor of luminosity distance fluctuations that was a linear 
combination of every eigenvector. 

Uncorrelated errors of apparent magnitudes were ran- 
domly chosen from a uniform distribution between 0.05 



and 0.2 mag. Moreover we assumed that any errors in the 
redshift are negligible. Uncorrelated magnitude and pe- 
culiar velocity scatter had fixed values of a m = 0.1 mag 
and a v = 300 km/s. 



B. Optimal quadratic estimator 

Because of the simplified form of our correlation matrix 
we chose the Newton iteration method for zero-finding of 
the derivative of probability function L, given by |27j 

lnL = - -^ln(dctC)-^pin2^, (20) 

where is A the data vector and C the correlation ma- 
trix given by Eq. ([9]). Our model depends on n p , p = 
1, . . . , N p parameters and is thus multidimensional. The 
correction to the parameters a p for Newton's method is 
then given by |2"Tll2l?] 



Sa v 



E 



d 2 lnL(q) 
ddpdcipi 



d\nL(a) 



da 



,pt 



(21) 



The common simplification is to replace the second 
derivative of In L with its ensemble average [571 HH] 

d 2 lni(a). 
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known as the Fisher matrix. When taking the ensemble 
average we assume that the underlying theory is correct 
and that the following relation is true (A A T ) = C. 
The correction to the parameter a p can then be written 

as [mug 

•fop - \ E( F_1 w Tr t( A AT - c ) ( c_1 c .p' c_1 )] 

(23) 

This method is called the optimal quadratic estimator 
(OQE) and is an iterative method. Although the OQE 
uses the Fisher matrix instead of the matrix of second 
derivatives, it converges to the same maximum. This is 
true because both matrices (F and the matrix of second 
derivatives) are invertible and for both we are in maxi- 
mum when 5a p = 0. The only approximation comes in 
using the Fisher matrix to approximate the errors. 

Since our parameters (P Q ,<T m , a v ) have physical mean- 
ing only when they are positive, we have checked this 
condition on every iteration step. Were they negative 
their values were put to zero. 



C. Application to synthetic data 

Equation Eq. ^ tells that aside from the matter power 
spectrum, the peculiar velocity correlation function de- 
pends on the cosmology through the linear growth factor 
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as well. Because we are interested in low-z SNe (z < 0.1), 
where the peculiar velocity effect is still large, we are in 
the regime where the Hubble rate is almost constant, and 
equal the Hubble constant. In the limit of low-z we can 
rewrite luminosity distances for flat universe as 
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With this in mind the factors in the luminosity distance 
correlation matrix (Eq. Q) become independent of cos- 
mological parameters. A little more care must be ex- 
ercised with the derivative of the growth factor. If we 
expand it into a more suitable form 
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Therefore, in this limit we can approximate 
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This illustrates that in the limit of small redshifts, the 
relevant quantity that our method is sensitive to is a 
non-dimensional growth factor multiplying matter power 
spectrum at some effective redshift (which we show to 
be the mean redshift). Our final result for the correla- 
tion matrix is thus algebraic combination of Eqs. (|8]), (|9|, 
m, pi p| and jm 



First, we have tested our method assuming a future 
survey of 10,000 supernovae la. We plot the results, 
together with the theoretical model used to create the 
dataset in Figure This illustrates that the method 
basically works. We will proceed with a series of tests. 

First, we have tested our method using synthetic data 
with different number of SN (1000,5000,10000) and differ- 
ent number of bins Nf, on which we estimated the matter 
power spectrum. We chose the bin positions by trial and 
error, so that the edge bin positions had very large errors 
and hence contained very little information. We therefore 
converged on the following set-up: one bin on small scales 
1 — 10/i/Mpc, two on large scales 10~ 4 — 10 _3 /i/Mpc and 
10~ 3 — 10~ 2 /i/Mpc. The rest were uniformly distributed 
in the logarithmic scale over the interval 0.01 — 1/i/Mpc. 
We always discard edge bins which contain essentially no 
information. 

Our method relies on the overall signal-to-noise to be 
high enough for the central theorem limit making the 
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FIG. 3: Shows the result of the new method in the 
form (f 2 P(k, z)) a (datapoints with errors) and the 
underlying model that went into creation of synthetic 
dataset (solid line). 



power spectrum constraints have Gaussian errors. To 
test this, we have performed a simple test. We analysed 
the same synthetic data two ways. First we measured 
the effective power spectrum from the data and then, for 
a standard AC DM cosmology measured the value of er 8 , 
while keeping all the other cosmological parameters fixed 
at their fiducial value. We then repeated the same, but 
this time skipped the intermediate test and measured a& 
directly from the synthetic data. In the limit of infinite 
number of supernovae, the two should match perfectly, 
but we do expect deviations to appear for a small finite 
number of datapoints. 

Results are shown in the Figure Q , where we plot the 
probability distributions for <jg for different number of 
SNe. On the vertical line we plotted the relative proba- 
bility L. As expected, with increasing number of SNe the 
distributions become more alike. If we fitted Gaussian 
distributions we found that the variances of the distri- 
butions for N sn = 5000 differ by roughly 40%, while the 
variances for N sn = 10000 differ only by a few percent 
(~ 4%). We stress that the point of this exercise is to 
see how many supernovae are required for reaching the 
Gaussian limit and that in general one should marginalise 
over other parameters. We also note that for any given 
realization, the maximum likelihood is expected to be 
distributed around the fiducial value according to the 
measurement error. 

Using the same method we checked how many bins are 
required. We plot this in the Figure (5). The distribu- 
tions vary in shape a little when changing number of bins 
as well as their positions in the k-space. But as long as 
we fix the number of bins around Nj, ~ 10 the variances 
of the distributions vary for only a few percent. 

We now turn to the dependence on the assumption of 
the fiducial cosmological model. To this end we generate 
synthetic data with cosmology A (Cl m = 0.24, w = — 1) 
and reconstruct the matter power spectrum P a with the 
same cosmology. Then we reconstruct the P a with differ- 
ent cosmology B (B ^ A) and compare the results with 
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FIG. 4: Show the probability distribution of 
cosmological parameter a& for data P a with respect to 

the reference power spectrum. It shows lines 
representing different number of SNe at the constant 
number of bins JVj = 9. We compare new method (red) 
with the brute- force calculation using Eq. ( 20 1 without 
approximations (black). Different line-styles represent 

different number of SNe: N sn = 1000 (solid), 
N sn = 5000 (dashed) and N sn = 10000 (dotted). P a 
stands for P a — (f 2 (z)P(k, z)) a . The normalization Lq 

is chosen such that L/Lq peaks at around 1. 



suit significantly and that was the purpose of this new 
method. In this way the most cosmological information 
is encoded in the result (((/ D) 2 P) a ). 

The last cross-check that we perform is to check if re- 
sult can indeed be approximated as if all supernovae were 
at a mean redshift. The matter power spectrum and / 
indeed vary with the redshift. In order to proceed, we di- 
vide the result of our method with average matter power 
spectrum for each bin (Pavg( z — z)), where z is the mean 
redshift of the input SNe synthetic data. If the quantity 
Pavg( z = z) is a good estimator of what we are measur- 
ing we should get a constant line at the value of f 2 . This 
is shown on figure ([7]), where we used 9 bins (JVj, = 9) 
and 10 realizations with N sn = 10, 000 SNe. We see that 
indeed the right / to use is the one evaluated at the mean 
redshift / = f(z — z). 

Figure Q shows that future surveys of several thou- 
sand SNe la will be able to estimate as to the accuracy 
of 6% using 10, 000 SNe la. 

Finally, we have investigated to what redshift we can 
push our method, by creating mock catalogs with increas- 
ing maximum redshfit. The approximations employed in 
this work clearly start to break down at z=0.2 and there- 
fore the data need to be split into several redshift bins if 
working over z = 0.1. 




FIG. 5: Shows the probability distribution of 
cosmological parameter <r 8 for data P a with respect to 
the reference power spectrum. Different line colours 
represent different number of bins at the constant 
number of SNe N sn = 10000: N b = 9 (red), N b = 5 
(blue) and N b — 12 (green). Black colour represents the 
brute-force calculation without using any 
approximations. P a stands for P a = (_/ (z)P(k,z)) a . 
The normalization Lq is chosen such that L/Lq peaks 
at around 1. 



respect to the error on the P a bins reconstructed with A. 

Figure [6] shows that the results vary with the input cos- 
mology only in the order of a few percent of the error and 
are thus negligible, when we vary parameter Q m and w 
within their present error-bars. Which shows that even if 
we use the slightly wrong fiducial cosmology in the power 
spectrum recovery procedure, this does not affect the re- 



V. CONCLUSIONS 

To sum up, we presented a new data compression 
method for cosmological extraction of data from the low- 
z SNe. We have shown that these supernovae effectively 
measure P(k, z)f(z) 2 at the mean redshift of the super- 
novae survey and that in the limit of large number of 
supernova N ~ 10, 000, the method is essentially op- 
timal and unbiased. We have further shown that the 
correlations of supernovae velocities are sourced by per- 
turbations on scales with wave-vectors ~ 10~ 2 Mpc//i to 
~ 0.5Mpc//i as these scales are the most constrained. 

Using mock catalogues we have shown that accuracy 
of 6% in as can be obtained using 10, 000 SNe la with 
the data compression method presented in this paper. 
We would also like to stress that using velocity angu- 
lar power spectrum should be excercised with caution, 
since CVs at different redshift bins are not correlated as 



was shown in section III Moreover, the data compression 
method presented in this paper does not carry any new 
information, it just shows that all information is repre- 
sented by f 2 (z) P(k, z). 

These basic conclusions are unlikely to change with 
the real data. However, at the limit of this number of su- 
pernovae, one should worry about potential systematics 
and biases. Two of those are the most important. First, 
the actual perturbations sourcing supernovae are not in 
the completely linear regime any more and second, su- 
pernovae don't trace velocity at random positions in the 
Universe, but instead trace it at position of rare density 



7 



— (n» 


,tu)s = (0.27, 


-1) 




,u>)s = (D.21, 











£ [h/Mpc] 



10" 



cq 
I 

ftf IB" 





(!]„ 


, ty, n t )s = (0-24, 


-1,0.01) "'"■■■:./ 




,w,fi*)B = (0-27. 


-1,0.01) 


■ (U, 


,a,n k ) B = (0.21, 


-1,0.01) 



fc [h/Mpc 




FIG. 6: Shows the difference of the P a bins at different cosmologies (A and B) with respect to the error on the 
result reconstructed with A. On the upper left plot we change Q m while keeping w constant (Q m = 0.27 - full line, 
fi OT = 0.21 - dashed; cosmology A had O m = 0.24). On the lower left plot we change w while keeping fi TO constant 

(w = —1.1 - full line, w = —0.9 - dashed; cosmology A had w = —1.0). The right-hand side plots are the same as 
left, but we have additionally used fi& = 0.01 in the assumed model. 
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FIG. 7: Shows the results of the new data compression 
method in the form (/ 2 P(fc, z)) a /P™f(z = z) on the k a 
bins. Two dashed constant lines represent values of / 
evaluated at two different redshifts [z = - green and 
z = z - blue; where z stands for the mean redshift of the 
synthetic SNe data). The mean of the results is shown 
as a full red line with ler errors (gray shaded region). 

The results are consistent with the latter 
(/ = f( z = z)). We used the OQE with N b = 9 bins 
and 10 realizations with N m = 10, 000 SNe. 



peaks in the primordial field, where galaxies eventually 
form. Corrections due to these effects can only be as- 
sessed by running iV-body simulations. This exceeds to 
scope of this paper. 

Recent SNe data contain of the order of a couple hun- 
dred low-redshift supernovae and so this method remains 
inapplicable. Unfortunately, the 10, 000 low-redshift SNe 
will not be available for the foreseeable future. The LSST 
survey predicts ~ 100, 000 SNe per year with redshift 
z < 1.2. Of this ~ 1400 SNe per year with z < 0.1. How- 
ever, there is no planned spectroscopic followup to these 
measurements. The LSST photometric redshift errors are 
predicted to be around 0.02(1 + z). We examined if in- 
formation on peculiar velocities could be extracted with 
such large redshift errors by approximating them with 
an extra velocity scatter of 6000km/s, but the results are 
not encouraging. Photometric errors of course have an 
additional issue of being highly non-Gaussian. 

If one could, however, measure 10,000 supernova la at 
low redshift, it might results in some interesting science. 
The auto-power spectrum of the population of galaxies 
in which supernova reside would measure the value of 
P(k)b 2 , where b is the galaxy bias, while the redshift- 
space distortions of the same sample would measure /3 = 
f/b. Combined with peculiar velocities measurements 
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discussed here, one would get an over-constrained system 
and if analysis was done on the same volume, the sample 
variance would cancel out. The quantity 



a 1 Pg 

S °-Wp7 n 



(28) 



would therefore have a unity value in Einstein gravity and 
deviation from unity would indicate either new physics or 
systematics. Note that in using several tracers over the 
same field, the sample variance would cancel and thus 
provide a much more stringent result than one might 
naively guess. Similar method using weak gravitational 
lensing was recently used in ([30 ). 



where k = |k|, x = |x| and 7 the angle between k in x. 
Also the following identity is true 



k d 
ix dk 



^ikx cos 7 ifcx cos 7 



fccos7. 



(A4) 



Now we are left with the integral of the exponent func- 
tion over the spherical harmonics. With the use of the 
following mathematical identities [3Tj 



^ikr cos 7 



^ a„j n (kx)P n {cosj), 



(A5) 



71=0 



ACKNOWLEDGEMENTS 



where a n = i n (2n + 1), j n are spherical Bessel functions 
and P n Legendre polynomials given by 



VI acknowledges support of the Berkeley Center of 
Cosmological Physics, where parts of this work were com- 
pleted during summer working visit. AS is supported in 
part by the U.S. Department of Energy under Contract 
No. DE-AC02-98CH10886. 



+n 



?fc, 9k 



(A6) 



where Y np are spherical harmonics, we can rewrite 



Eq. (A3 I into 



APPENDIX 



When computing the angular power spectrum we get 
(ae m a* e , m ,) = D'(z 1 )D'(z 2 ) 



^p(k) f dn k - 



fc 4 



(27T) 

I ^i)^/ m (iiK k xi k-x r 

J ^(x 2 )^ m ,(x 2 ) e - lk X2 k-x 2 , (Al) 



where ag m stand for the expansion coefficients, D(z) is 
normalized linear growth factor and derivatives are with 
respect to conformal time {a{t)drj = cdt). In equation 
(Al) there are two integrals of the form 



i tm = J rfO(x)r; m (x)e Ikx k-x, 



(A2) 



where we have dropped the SN label in the subs cript 
(1 or 2). We notice that the integrals in Eq. (Al) 
are complex conjugates, so we only need to calcu- 
late one. If we choose the coordinate system such 
that k = k (sin 9 k cos <pk, sin 6k sinc/^, cos Ok) and x = 
x (sin 6 cos (f>, sin sin 0, cos 9) we can rewrite the above 
integral into 



i tm = / dfty/ m (MK fexcoS7 fccos 7 , 



(A3) 



4"7T 

CL-n. ~ ~ jn ( ) ' 



m ix dk J ~" 2n + 1 

n—0 

+ 11 



Yn P * (9 k , M J dQY; m (e^)Y np (9, 



p——n 



(A7) 



where the only integral left is the orthonormal relation 
of spherical harmonics 



/ 



dftY£ m (9, 4>)Y np (9, 4>) - u tn u mp , 



(A8) 



where 5's on the right-hand side are Kronecker delta's. 
Using the above identity we can evaluate both sums in 
the Eq. (|A7|) and write the result as 



k d 4-7T 
hm = — -Er.ae —„ ■ ji(kx)Y^ m (0 k ,4> k ) 



ixdk 2£+V 
4?rfcl xdk Yh 



(A9) 



Comparing the result with the original equation ( Al |, we 



have computed Ii m (kxi). With complex conjugation and 
variable substitution of this result we can write the sec- 
ond integral as Ip m ,(kx2), where x± t 2 — | x i,2|- Rewriting 



9 



Eq. ( Al I and inserting the integral values we get which is exactly what we were looking for. The term on 

the right-hand side of the equation, that stands in front 



{ai m a e ,. 



D'(z )D'(z ) f fc2rffc P(k) ° f the Su ' Smm ' is by dennition equal to C t . Implying the 

J (2w) 3 Kronecker delta 8u> we can write 



dQ k -^Ii m (kx 1 )I e/m ,(kx 2 ) 
=D'( Zl )D'(z 2 ) f ^P(fc)- 



du{kx 1 )dj tl {kx 2 ) . l _ l ' C l {z 1 ,z 2 )=D'(z 1 )D'{z 2 ) 
x\dk x 2 dk 



d& k Y? m (8 k ,4> k )Y Vm ,(6 h ,<l> k ) ■ I dkP ( k ) 



kie m \ c 'k,<Pk)ie> m '{Ok,<Pk) J Q w y Xl Qk J\ x 2 dk 



'-D\ Zl )D'(z 2 ) J ^P(fc)- 
dji(k Xl ) dj t (kx 2 ) . e _ e > 

' ?T, 7^, 1 OU'Omm'i 

X\Ok x 2 ok 

(A10) 



2 

7T 

djt(kxi)\ ( dj e (kx 2 ) 



(All) 
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